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ABSTRACT 

Multifractal scaling (MFS) refers to structures that can be described as a collection 
of interwoven fractal subsets which exhibit power-law spatial scaling behavior with a 
range of scaling exponents (concentration, or singularity, strengths) and dimensions. 
The existence of MFS implies an underlying multiplicative (or hierarchical, or cascade) 
process. Panoramic column density images of several nearby star-forming cloud 
complexes, constructed from IRAS data, are shown to exhibit such multifractal 
scaling, which we interpret as indirect but quantitative evidence for nested hierarchical 
structure. The relation between the dimensions of the subsets and their concentration 
strengths (the "multifractal spectrum") appears to satisfactorily order the observed 
regions in terms of the mixture of geometries present, from strong point-like 
concentrations, to line-like filaments or fronts, to space-filling diffuse structures. This 
multifractal spectrum is a global property of the regions studied, and does not rely on 
any operational definition of "clouds." The range of forms of the multifractal spectrum 
among the regions studied implies that the column density structures do not form 
a universality class, in contrast to indications for velocity and passive scalar fields 
in incompressible turbulence, providing another indication that the physics of highly 
compressible interstellar gas dynamics differs fundamentally from incompressible 
turbulence. There is no correlation between the geometrical properties of the regions 
studied and their level of internal star formation activity, a result that is also apparent 
from visual inspection. We discuss the viability of the multifractal spectrum as a 
measure of the structural "complexity" of the regions studied, and emphasize the 
problematic dependence of all structural descriptors on the subjective pre-selection 
of the region to be described. A comparison of IRAS 100 fim column density (not 
intensity) images with ^^CO, visual extinction, and C^^O data suggests that structural 
details are captured by IRAS up to at least 30 magnitudes of visual extinction, except 
in the vicinity of embedded stars, and that lower-column density connective structure 
not seen by other methods is revealed. 



Subject headings: galaxies: ISM — ISM: clouds, structure — methods: data analysis 
— techniques: image processing — turbulence 



1. Introduction 

Multifractal scaling (MFS) is a property that refers to measures or fields for which the local 
intensity (e.g. density, column density, radial velocity, temperature), integrated out to distance r 
from a given position, scales as a power law r", with the scaling exponent a varying from position 
to position. If the collection of positions with scaling exponents in (a, a + da) can be described as 
a fractal with (Hausdorff ) dimension /(a) , for a range of a, then the structure is multifractal, and 
the function f(a) is called the multifractal spectrum (MFS; we use MFS for both "multifractal 
scaling" and "multifractal spectrum" when no ambiguity can arise). This definition is not unique 
in the literature, but is sufficient for our purposes (see Feder 1988, McCawley 1990, Muzzio 
et al. 1992). The exponent a is called the "singularity strength" in the multifractal literature, since 
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the power law in principle gives a density singularity as r approaches zero; here we prefer to use 
the more descriptive term "concentration strength" for this exponent. Similarly, the multifractal 
spectrum is often referred to as the "spectrum of singularities" in the physics hterature, but we 
avoid this terminology here; however, since we do use the potentially misleading term "multifractal 
spectrum" , it is important to remember that "/" is a dimension of a subset, not a filling factor. 

The MFS basically contains information on how the geometry, as measured by the dimension, 
depends on the concentration strength a; loosely speaking, the /(a) curve shows the manner 
in which regions of different concentration strength fill space. For example, for subregions with 
a given a, the dimension will be different for point-like, line (or filament) -like, or space-filling 
structures. So the MFS indicates the distribution of geometries present in a complex structure. 
For a homogeneous fractal the /(a) curve collapses to a point; otherwise the measure (intensity 
image) is multifractal, since regions of different concentration strength have different dimensions. 
A common /(o) curve between physical processes and/or mathematical models can be interpreted 
as evidence that they belong to the same "universality class" ,[| at least with respect to the 
f{d) function. However, it is well-known that the /(a) descriptor (and all other structural 
descriptors of which we are aware) are degenerate with respect to the underlying physics, in the 
sense that systems with very different underlying dynamics can exhibit identical /(o) curves within 
the measurement uncertainties; see Bulajich and Perez-Pascual (1991), also Chhabra, Jensen, and 
Sreenivasan (1989). We give a simple example in Appendix B. 

Multifractal scaling usually occurs in systems that arise from multiplicative processes, or 
"cascades," in which the intensity (e.g. density, etc.) at a given position is due to successive 
multiplications along a hierarchical tree (see below), so the existence of MFS (and its properties) 
can be used to infer hierarchical (nested) structure, and to constrain the types of physical processes 
that give rise to such a structure. The present work attempts to test for MFS in column density 
maps of nearby cloud complexes and to find whether the resulting f{d) function can be used to 
describe the mixture of geometries present in these structures. 

Densely sampled extinction, molecular line, H I line, and IRAS maps of nearby cloud 
complexes have shown that the gas and dust are organized into complicated structures, with 
irregularities present over a wide range of scales, from roughly 0.02 parsecs (resolution limit for the 
nearest cloud complexes) to at least a few hundred parsecs (see reviews in Scalo 1985, 1990, Wilson 
and Walmsley 1989, Falgarone and Phillips 1991, Stutzki et al. 1991, Genzel 1991, Falgarone and 
Perault 1992, Elmegreen 1992, and references therein). The geometry and scaling of this structure 
may provide a signature of the physical processes that are involved in star formation, and may 
be related to important quantities such as the star formation efficiency and the initial mass 
function, as suggested by Henriksen (1991) and Larson (1992). Elmegreen and Falgarone (1996) 
and Elmegreen (1997) have proposed that the fractal scaling properties of the interstellar medium 
(ISM) are consistent with many of its observed features, including the mass-size scaling, the cloud 
mass spectrum, and the existence and properties of the intercloud medium. Scaling behavior is 
the basis for physical interpretation in a huge number of applications within physics, with phase 
transitions and turbulence being most prominent. 

Since there is no unique way to characterize structure, however, there are many ways to 
define structural scaling. One such method is the perimeter-area relation (see Feder 1988, chap. 
12). If contours of equal intensity (or column density, or temperature, etc.) exhibit a power law 
perimeter-area relation with a non-integer exponent over a range of scales, this exponent may 
be interpreted in terms of a fractal dimension which characterizes the manner in which these 
curves fill space. A number of studies have demonstrated such behavior in local interstellar cloud 



^ For example, in turbulence the MFS is believed to be "universal" at the onset of chaos, because it agrees with 
the MFS of the 1-dimensional circle map. 



- 3 - 



structures structures^ using extinction maps (Beech 1987, Hetem and Lepine 1993), H I maps 
(Wakker 1990), 100 /xm dust intensity or column density maps (Bazell and Desert 1988, Scalo 1990, 
Dickman, Horvath, and Margulis 1990, Vogelaar, Wakker, and Schwarz 1991) and CO emission 
maps (Falgarone, Phillips, and Walker 1991); see also Vogelaar and Wakker (1994). Surprisingly, 
these studies find a perimeter-area dimension of about 1.3-1.5 for most clouds studied, similar 
to the value found for terrestrial cloud and rain areas (Lovejoy 1982) and slices of laboratory 
turbulence (see Sreenivasan 1991 for a review). Even the HI gas in the M81 group of galaxies 
appears to exhibit a very similar perimeter-area dimension (Westpfahl et al 1999). This result 
suggests a sort of universality for cloud fractal geometry which is independent of the presence of 
young stars or the importance of self-gravity, and may be related to turbulence (by which we mean 
phenomena associated with nonlinear advection), as emphasized by Scalo (1990) and Falgarone 
et al. (1991), or some other process, since the dimension of a fractal is extremely degenerate with 
respect to underlying processes and morphological appearance (i.e. many very different processes 
and forms can give rise to the same dimension). 

In the present work we are concerned with a different type of scaling behavior, and it is 
important to appreciate the difference. If the intensity (or column density, or radial velocity, 
or any measured quantity) image is represented as a surface whose local height is the intensity 
("mountain range," or "topographic," representation), then the dimension obtained from the 
perimeter-area relation characterizes the irregularity of contours formed by horizontal cuts through 
the "mountain range." However, this dimension does not uniquely or completely characterize 
the scaling of the intensity surface itself. Adding a constant or multiplying by any function of 
the image intensity will not affect the power-law slope of the perimeter-area relation, although, 
for example, the "spikiness" of the mountain range may be greatly affected. Furthermore, it 
is commonly stated that if the perimeter-area relation is a power law with exponent d, then 
the clouds have a fractal dimension D = 2d. However, this gives the incorrect impression that 
if this scaling holds then the cloud complex is a statistically self-similar fractal, when in fact 

perimeter-area scaling only demonstrates that the "mountain range" is a self-affine fractal.^ The 
structure of the intensity surface of such self-afRne fractals may possess a continuous range of 
scaling exponents or dimensions. 

In the multifractal approach used here, the column density "surface" is not sliced at different 
intensities. Instead the present technique is equivalent to partitioning the image into different 
subsets according to the value of the local concentration strength a, defined by the exponent of 
a local power law scaling of the intensity integrated out to distance r from a position, if such 
a power-law scaling exists. The dimension of this subset is not estimated by treating it as a 
collection of curves (contour map) and characterizing the irregularity of each curve; the dimension 
is estimated for the entire subset, treated as a collection of points, and includes all the points, not 
just those comprising disjoint closed curves. No reduction to "contours" is involved in the process. 
As will be seen, the dimensions estimated in this way are very different from the perimeter-area 
dimension, although the relation between the two remains unclear to us. For example, we will 
show that the column density surface in the Taurus complex (and the other regions we examined) 
is a multifractal, even though it is known to possess a well-defined perimeter-area relation with 
a single fractal dimension for all intensities (Scalo 1990, Falgarone et al. 1991). Still, we think 
that the perimeter-area dimension is a useful and illuminating tool that is complementary to 
other descriptors. A detailed discussion of the perimeter-area dimension applied to the interstellar 
medium is given in Dickman et al. (1990), to which the reader is referred for more discussion. 



^ We use the term "cloud" loosely to designate an under-resolved observed region of space that appears to contain 
a coherent density enhancement, recognizing that under sufficient resolution this apparent entity will likely lose its 
coherence and break up into complex substructure. 

Self-similarity refers to structures for which each sub-structure can be obtained from the whole structure by a 
linear contraction which reduces the original structure by the same scale factor in all coordinates. Self-afhne refers 
to structures whose substructure can be obtained from the original by a linear transformation in which different 
coordinates (e.g. column density and spatial scale) are contracted by different factors. 
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(Note, however, that the latter study was based on IRAS intensity images, not column densities 
as used in the present MFS investigation.) 

It is also important to distinguish the structural scaling examined here from the type of power 
law correlations which have often been discussed for interstellar clouds, e.g. cloud size or mass 
spectrum, or correlations between density or velocity dispersion and region size. First, these latter 
correlations are between properties defined for separate entities, such as "clumps within clouds" or 
certain operationally-defined types of "clouds" located in different star-forming regions or different 
locations in the galaxy. In contrast, the scaling examined here is a property of a given region as 
a whole, not a correlation between attributes of any operational entities ("clouds"); the concept 
of "cloud" does not enter the scaling examined here. Second, the specific form of the quantities 
whose scaling is being examined here (see Section III below) is directly related to the geometry of 
the structure, in the sense that they describe how the structure fills space. 

The manner in which the discovery of multifractal scaling might be related to physical 
processes is indicated by similar investigations in other fields. Fully-developed incompressible 
turbulence in a variety of experimental situations exhibits an apparently universal multifractal 
singularity spectrum, whether determined for the spatial distribution of the dissipation field 
(see Sreenivasan 1991, 1996 and references therein) or for the velocity signal (Muzy, Bacry, and 
Arneodo 1991). Thus the /(a) spectrum may reveal an invariant signature of the underlying 
nonlinear fiuid dynamical processes in incompressible turbulence (e.g. vorticity tilting and 
stretching in three dimensions); i.e. f{d) may reflect the essential physical processes independently 
of the type of flow or boundary conditions. However, the process reflected by the MFS may be 
much more generic. For example, Meneveau and Sreenivasan (1989), Chhabra and Sreenivasan 
(1992) and others have argued that turbulence is in the same MFS universality class as certain 
simple multiplicative processes, providing a guide for theoretical modeling. A comparison with 
direct numerical simulations was given in Hosokawa and Yamamoto (1990). Here we show that 
such a universal /(a) function does not exist for the density structures of the highly compressible 
cool interstellar medium, even though the structures are multifractal. 

In an extragalactic context, Martinez et al. (1990) suggested that the three-dimensional 
distribution of galaxies in the CfA survey possesses multifractal scaling, and showed that the 
corresponding f{a) spectrum could be used to distinguish the Soneira- Peebles hierarchical model 
(a homogeneous fractal) from the Voronoi tessellation model (a multifractal), even though the 
correlation function for these two models are very similar. A study of multifractal scaling in 
cold dark matter simulations, with an emphasis on finite size effects, has been given by Colombi 
et al. (1992). A number of more recent discussions of the relation between possible multifractal 
scaling and large-scale galaxy clustering have appeared; see Sylos Labini and Pietronero (1996) 
and references given there. Multifractal scaling has been discussed in other astrophysical contexts, 
such as the attempt to relate the multifractal spectrum to the geometry and underlying physics 
of the solar photospheric magnetic field by Lawrence, Rumaikin, and Cadavid (1993). Outside of 
astrophysics, the use of the multifractal spectrum as a diagnostic tool is widespread. Examples 
include various aggregation and coalescence models (see Nagatani 1992), characterization of stages 
of cancer development through analysis of normal and malignant tissues (Muller et al. 1991), and 
analysis of protein folding through profiles of solvent accessibilities of amino acid chains (Balafas 
and Dewey 1995). 

Multifractal scaling is closely associated with multiplicative processes (see sec. II below), 
such as stretching and folding of fluid elements (e.g. Muzzio et al. 1992), and so its existence is 
sometimes taken as evidence for hierarchical structure (e.g. Balafas and Dewey 1995). Since the 
hierarchical nature of interstellar cloud structure is not yet adequately established or characterized 



In contrast, the existence of other types of power-law scahng does not obviously provide evidence for nested or 
multiplicative structure. For example, if a collection of disks with a power law distribution of sizes is distributed at 
random on a plane, the distribution is a homogeneous fractal, even though there is no nesting; i.e. the distribution 
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quantitatively (e.g. Houlahan and Scalo 1992), a search for multifractal scaling in ISM data is of 
obvious interest. In particular, Meneveau and Chhabra (1990) and O'Brian and Chhabra (1992) 
have suggested an empirical test for hierarchical processes based on the MFS, which we apply 
below. 

In addition, the MFS may be useful in characterizing complex interstellar structures in terms 
of the distribution of geometries, since it describes the dimensionality of the subsets as a function 
of concentration strength. As far as we know, other global descriptors of structure (e.g. the 
correlation function) are insensitive to geometry, except perhaps for fairly simple structures, like 
periodically-spaced density fluctuations along a filament. Descriptors can be found to identify 
local specific forms (e.g. minimal spanning tree for finding filaments, Bhavsar and Ling 1988; 
line and plane shape detectors based on moments, Luo, Vishniac, and Martel 1996), but these 
descriptors do not characterize the global geometry in the manner of the MFS. 



2. Data 

We obtained 60 and 100 fj,m IRAS images of four different regions of low-mass star formation 
from the Infrared Processing and Analysis Center (IPAC). The emission seen at these wavelengths 

is the radiation of warm dust grains heated by local stars and the galactic radiation field. IPAC 
estimated and subtracted the local zodiacal background, and we removed the galactic background 
using a cosecant model profile. Other background models, such as fitting a plane or polynomial to 
local minima, yielded very similar results. A detailed discussion of these and other points can be 
found in Verter and Rickard (1998). 

The far-infrared emission depends on both the dust temperature and the column density. 
These two quantities may be separated given the 60 and 100 micron fluxes if it is assumed 
that all of the dust along a given line of sight is at the same temperature, the dust is optically 
thin, and that a power-law dust emissivity holds. In such a model, the infrared flux is given by 
Fx = Bx{Taust){^ — e~^) where r is the dust optical depth and BxiT^ust) is the Planck function. 
The dust optical depth was assumed to vary with wavelength as A", with n = -1 or -2. This 
relation may be iteratively solved for the temperature and column density for every point in the 
image. 

This method for estimating dust temperatures and column densities assumes that the 
dust along a given line of sight is at a constant temperature, which biases the column density 
estimates toward the column density of the warmer dust. For embedded sources this warm dust 
will be concentrated very near the source (compared to the line of sight depth), constituting a 
small fraction of the true column density along the line of sight, so the column density will be 
underestimated, sometimes by a large factor. As a result, embedded protostars and H II regions 
may appear as holes in the final images. (These artifactual holes may actually provide a useful 
method for identifying embedded sources, since they are visually very distinct.) For this reason we 
mostly chose regions of relatively low mass star formation for this study, since there the stars have 
minimal effect on the column density images. A more sophisticated technique that provides an 
estimate of temperature variations along the line of sight has been proposed by Xie et al. (1991). 

The resulting maps do, however, strongly resemble column density maps based on extinction 
and molecular line data away from embedded stars (e.g. Langer et al. 1989 for the B5 cloud, 
Snell et al. 1989 for Heiles Cloud 2 and B18 in Taurus, Jarrett, Dickman, and Herbst 1990 for p 
Oph; below and Appendix A for various subregions in Taurus) . We emphasize that this agreement 
applies to column density maps, and that the intensity maps can appear significantly different; 
some previous comparisons of IRAS maps with molecular line (Abergel et al. 1995) or extinction 



will appear (statistically) self-similar as one "zooms in" on it, but the structure is not hierarchical with respect to 
nesting. There are many ways to obtain a power law size distribution that do not involve any multiplicative processes. 
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(Arcc & Goodman 1999) images have used 60 and 100 fim intensity images. The rehabihty of 
the column density maps presented here can best be seen by visual comparison with the recent 
adaptive grid wavelet transform extinction maps for the same regions given by Cambresy (1999). 

Four nearby (~ 100-200 pc) cloud complexes containing low- to intermediate-mass young 

stellar objects were chosen for study: Taurus, R Corona Australis, Chamaeleon-Musca, and 
Scorpius-Ophiuchus. The last three regions were divided into 2, 4 and 6 subregions on the basis 
of their differing visual appearance and star formation properties. Table 1 lists the positions, 
estimated distances, number of resolution elements (assuming a resolution of 3 arcmin), and rough 
estimates of the projected density of young stellar objects, when available, in each subregion. 
Molecular line maps of these regions with the NANTEN telescope are summarized in Kawamura 
et al. (1999) and Yonekura ct al. (2000); comparisons with earlier work and additional references 
can be found there. Extinction maps that cover roughly the same areas as examined here can 
be found in Cernicharo (1991, Taurus), Rossano (1978a, R Cr A), Gregorio-Hetem et al. (1988, 
Cham-Musca), Rossano (1978b, Sco-Oph), Arce & Goodman (1999, Taurus), and especially 
Cambresy (1999, Taurus, Sco-Oph, Cham-Musca, and R Cr A). The star formation activity 
varies from no detected activity in Cham III to widely distributed low-mass stars (and probably 
higher-mass B stars — see Walter and Boyd 1991) in Taurus, to dense and probably bound clusters 
in the p Oph core and in the R Cr A core. A discussion of the young stellar populations in parts of 
these regions with references to the original papers can be found in the review by Zinnecker et al. 
(1993) and the study by Chen et al. (1997). The gas in the entire Sco-Oph region has been affected 
by interaction with massive stars (see DeGcus 1988). Figures 1, 2, and 3 show logarithmically 
scaled column density maps of the Cham-Musca, RCrA, and Sco-Oph regions. An image of Taurus 
enhanced by unsharp masking is shown in Appendix A. The original Taurus column density image 
is given in Scalo (1990) and a contour map and a "landscape" representation of the region are in 
Houlahan and Scalo (1992). 




Fig. 1. — Logarithmically scaled image of the column density distribution Chamaeleon-Musca. The 
small dark holes arc embedded protostars which locally heat the dust, yielding an underestimate 
of the actual column density. The image is a mosaic of three frames, each of which covers 8 x 8° on 
the sky, corresponding to 20 x 20pc (assuming a distance of 140 pc), and contains roughly 6.7 x 10^ 
resolution elements. North is up and west is to the right. 
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Fig. 2. — Logarithmically scaled image of the column density distribution around R Corona 

Australis. The image covers 14 x 8° on the sky corresponding to 32x18 pc (assuming a distance of 
130 pc) and contains roughly 4.5x10^ resolution elements. Residual striping artifacts run vertically 
through the R CrA W subregion. North is up and west is to the right. 




Fig. 3. — Logarithmically scaled image of the column density distribution in the Scorpius-Ophiuchus 
region. The image covers roughly 28 x 28° on the sky corresponding to 61x61 pc (assuming a 
distance of 125 pc), and contains 2.5x10^ resolution elements. The light and dark bands running 
below and to the right of L134 and east of p Oph are residual striping artifacts. North is up and 
west is to the right. 

There is considerable uncertainty concerning the degree to which IRAS 100 /xm column 
density Nioo represents the true column density. A thorough discussion, with application to cirrus 
clouds, has been given by Verter and Rickard (1998). Generally the relation between Nioo and 
^^CO column density exhibits a great deal of scatter, yet the resulting spatal images of Nioo do 
closely resemble the structure of ^^CO and extinction maps (e.g. Langer et al. 1989 for the B5 
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cloud, Snell et al. 1989 for Heiles Cloud 2 and B18 in Taurus, Jarret, Dickman, and Hcrbst 1990 
for p Oph). Much of the scatter probably arises from a combination of the residual striping, the 
embedded protostar effect discussed above, and perhaps the background subtraction at the lowest 
intensity levels. (Wc point out that the large scatter found in studies like Abcrgel et al. 1995 is 
most likely due to the use of IRAS intensities rather than column densities.) However none of the 
effects seriously compromises the overall spatial structure except in very localized ways, or at the 
lowest intensity levels. The latter effect forces us to consider only positive moments in the MFS 
(i.e. the "hole" geometry is not considered - see sec. 4.1 below). 

Since the MFS responds to the global structure on all scales, it is the similarity between the 
spatial structure in Nioo and other tracers that is most important. Wc have compared the images 
of Nioo with available ^^CO and extinction maps of five subregions in Taurus, and the results 
show excellent agreement. An example, the subregion L1529 filamentary substructure in Taurus is 
shown in Appendix A. The only serious disagreement occurs for the very opaque L1495 subregion. 
However we find that the Nioo image for this subregion docs agree much better with the C^^O 
image of Duvert, Ccrnicharo, and Baudry (1986). This suggests that Nioo is a better tracer of 
column density than either ^^CO or visual extinction, both of which have severe optical depth 
effects that set in around Ay ~ 3 — 5. The case of L1495 suggests that the upper limit on Ay for 
the validity of Nioo may be around 30 mag, rather than the 10 mag limit suggested by Jarrctt 
et al. (1989) for the 60 //m column density. We also find surprisingly good agreement between the 
p Oph core Nioo image derived here and the submillimeter image given by Motte et al. (1998), 
considering the much better resolution of the submillimeter image, suggesting that the Ay limit 
for Nioo (except along lines of sight to embedded protostars) may exceed 50 mag. Thus the IRAS 
column density upper limit may be similar to that of infrared extinction (e.g. Lada, Alves, and 
Lada 1999), except that IRAS is locally thwarted by embedded stars, while infrared extinction 
requires smoothing to reduce noise introduced by structure in the background star positions. 

A strong indication of the validity of the IRAS column density images is the excellent visual 
agreement of our images with the adaptive grid/wavelet decomposition optical extinction maps of 
several of the regions by Cambresy (1999). Another favorable example is that the "cores" mapped 
in the C^^O study of a large part of Ophiuchus by Tachihara, Mizuno, & Fukui (1999) can be seen 
to correspond to the highest-column density peaks in the IRAS images of Ophiuchus displayed in 
Fig.3; the IRAS column density images indicate that the apparently isolated C^^O cores studied 
by Tachihara et al., especially in the regions we refer to as Oph N and U, are actually connected in 
a web of elongated or filamentary forms whose column density was too small for the C^^O survey. 

Overall, these comparisons lead us to be optimistic about the validity of the IRAS column 
density maps, except near embedded sources or at very low column densities (which we do not use 
in the subsequent analysis). As remarked above, we found our results to be robust with respect to 
the adopted power law form of the far-infrared emissivity. Lastly, we found no problems associated 
with producing seamless mosaics of the IRAS maps, except that it was technically laborious. 

We emphasize that a major advantage of the IRAS data is that it can probe a very large 
range of column densities compared to other methods, roughly from about 2 x lO^^cm"^ (see 
Weiland eta 1 1986, Verter and Rickard 1998) to 6 x lO^^cm"^ (assuming that the ratio of Ay to 
rioo is a few thousand). 



3. Calculation of Multifractal Statistics 

The multifractal spectrum is a method of characterizing images and point sets which possess 
a type of statistical scale invariance called multifractal scaling. The local density profile around 
a point in the image is characterized by a scaling exponent called the singularity strength (here 
we use the term "concentration strength" ) . If Pj (L) is the integral of the image intensity (column 
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density in the present work) within an aperture of radius L centered at position i, then the 
concentration strength a is defined by Pi{L) oc L". For a two-dimensional homogeneous image 
a = 2, while for a sharp spike a = 0. Thus, the concentration strength measures the local 
"spikincss" in the image. The function /(a) measures how all the points in the image which have 
a given concentration strength are distributed and it is the Hausdorff dimension of that point set. 
The numerical calculation of /(a) has proven to be delicate. Difficulties with procedures used in 
multifractal analysis have been discussed by Chhabra, Jensen, & Sreenivasan (1989), Yamaguti 
and Prado (1995) and Veneziano, Moglin, and Bras (1995). The latter paper, in particular, 
emphasizes that the /(a) spectrum only generates the upper envelope of the true MFS, while 
missing interior points not in the upper envelope, and generates spurious points if the actual 
MFS is discrete (a bifractal would be an extreme example). We selected the canonical method 
of Chhabra and Jensen (1989), which is computationally efficient, avoids some of the problems 
faced by other techniques, and provides an intriguing mathematical connection to methods used 
in statistical mechanics. (See McCauley 1990 for a discussion of "thermodynamic formalism" for 
multifractals.) 

The canonical method for computing /(o) begins by partitioning the image into boxes of 
width L. If Pi{L) is the intensity integrated over box i, then /(a) is given implicitly by 

/(,) = hm ^^^-yi°g^-(g) (1) 
•' L->o log L ^ ' 

= = ^— (2) 

i ° 

where iJLi{q) = Pf{L)/Y^Pf are the q^^ coarse-grained moments of the original image and the 

i 

sums cover the entire image. For real data sets, the limits are estimated by least square fitting 
of J2fJ'i^og Hi and ^^ilogPi as a function of logL for different q values. If the linear fits are 

i i 

good and the slope varies with q, then the data set is said to possess multifractal scaling over the 
specified range of scales. If / and a take on only a single value for all moments q, then the image 
is a homogeneous fractal. 

The order parameter q can be shown to be the slope of the /(o) curve, so the /(a) curve for 
multifractals is a convex downward opening curve peaking at the point where q = 0. The value 
of / at that point is the Hausdorff dimension of the support of the measure, which is 2 for the 
images considered here, while a{q = 0) is just the average concentration strength of the image. At 
the point where q = 1, the f{a) curve must lie on the diagonal line /(a) =a. The point at g = cx) 
characterizes the regions in the image which have the smallest values of a which may be thought 
of as the sharpest spikes. Similarly, the q = — oo point characterizes the regions which have the 
steepest holes, i.e. with the largest values of a. The /(«) function is related to the so-called 
generalized or Renyi dimensions Dq (which characterize the scaling of the qth order partition 
function ^Pf{L)) through a Legendre transform. 

i 

The dependence of the MFS on the parameters of a simple "toy model" multiplicative process 
is illustrated in Appendix B. 

The referee has pointed out that the formalism involved in obtaining the MFS seems rather 
abstruse and non-intuitive, and suggests the general principle, with which we agree, that any 
empirical descriptor should bear a proportional requirement of success in terms of physical 
interpretation in exchange for complexity of effort. While we agree that the canonical formalism 
adopted here is not at all transparent, there are several points which suggest that the use of the 
MFS does not violate the above general principle. First, the computational complexity, if measured 
in cpu time, is actually smaller for the MFS than for the correlation function. Second, it must be 
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recognized that the degree to which a descriptor seems obscure or complex depends inversely on 
the familiarity of the user. For example, the meaning and significance of the correlation function 
often seems obscure to those who have not dealt directly with it. Concerning the ability of the 
MFS to reveal some underlying physics, the reader may judge from the rest of the paper. We only 
remark that the MFS does allow us to infer the presence of some multiplicative process, and also 
shows a fundamental difference between incompressible and supersonic turbulence with respect to 
universality class, as discussed below. By comparison, studies of the velocity correlation function, 
for example, have revealed few significant clues to underlying physical processes, as far as we know. 



4. Results 
4.1. Multifractal Scaling 

Figure 4 shows J2l^i logMi ^-^d X^/^i log-Pi as a function of log L for various positive q values in 

i i 

the Taurus region. The range of box sizes used for the fits was 10 to 100 pixels, which corresponds 
to linear sizes of 0.4 to 4 pc. The lower limit is set by the resolution. For example, the adopted 
IRAS angular resolution of 3 arcmin corresponds to a linear resolution of 0.13 pc for Taurus, 
at the assumed distance of 140 pc (similar to all the regions studied here). Image structure is 
generally severely compromised for images smaller than 3-4 resolution elements, or 0.3-0.5 pc for 
the regions studied here. The upper limit is set by the fact that at larger scales a significant 
fraction of the boxes used to compute the MFS overlap the edges of the images. We emphasize 
that the image size must he significantly larger than the upper limit for the range examined for 
scaling behavior; otherwise edge effects distort the scaling. Related edge effects give comparable 
distortion to other descriptors of structure (e.g. correlation function, see Houlahan and Scalo 
1990 for an explicit demonstration). The panoramic nature of the IRAS images make it possible 
to satisfy this requirement. However, errors in calibration and background removal dominate and 
distort the scaling for negative q values, which accentuate the low-flux areas. As a result, the 
scaling for negative q values could not be investigated. The good linear fits (at least for lower q 
values), with different slopes for different q values, exhibited in Fig. 4 indicate that the column 
density structure in Taurus is multifractal, at least over the range of scales examined. These fits 
are typical of nearly all the subregions. Fig. 5 shows the scaling behavior of all of the subregions 
for a q value of 8. Wc conclude that all the regions used in this study exhibit multifractal scaling. 
This result demonstrates that there are no preferred scales in the column density structure over a 
factor of about ten in scale (roughly 0.4 to 4 pc in the regions studied). Such a small range does 
not provide a very strong constraint on preferred scales, however. 

Multifractal scaling is usually associated with multiplicative processes, involving hierarchical 
nesting or a "cascade." Meneveau and Chhabra (1990) and O'Brian and Chhabra (1992) have 
shown that the two-point correlation function of the spatial field of scaling exponents a, Cq,(s), 
should exhibit a logarithmic decline with separation s for random isotropic multiplicative processes, 
providing a test for hierarchical structure. Wc have applied this test to the Taurus region by 
constructing an images of a, and then calculating the correlation function C„(s). The result is 
shown in Fig. 6, where it is seen that a significant range of logarithmic decay of Ca{s) exists. 
The fit to a logarithmic decay is even better than the turbulence examples shown in O'Brian and 
Chhabra (1992). We postpone a discussion of the other regions to a separate publication. 







I I I I I " 



1.2 1.4 1.6 1.8 
log(box size) 



1.2 1.4 1.6 1.1 
log(box size) 



Fig. 4. — Multifractal scaling in Taurus. The slopes of the linear regression fits (indicated by 
straight lines) in the left panel gives the value of a for each q value while the slopes in the right 
panel give the valves of /. The units of the abscissas are log arcminutes. At an assumed distance 
of 140 pc, the box size limits range from 0.4 to 4.0 pc. Arbitrary constants were added to each set 
of points to prevent the fits from overlapping. 




Fig. 5. — Multifractal scaling for all the regions. The fits for a and / at g = 8 are shown for all 
the regions. The fits are ordered by their slopes, with the steepest slopes at the top. Each set of 
data points were arbitrarily shifted up or down to prevent overlapping on the plot. 
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Fig. 6. — Correlation function for the image of concentration or "singularity" strengths a for the 
Taurus region. The significant range over which logarithmic decay occurs is interpreted as evidence 
for nested hierarchical structure (see text). 



4.2. Multifractal Spectrum 

Figure 7 shows the /(a) curve for positive q values for each of the subregions. These 
/(a) curves were found to be fairly insensitive to the choice of background or to the adopted power 
law in the variation of dust emissivity with wavelength. The error bars represent the fitting 
error in determining / and a. The curves are ordered by their width, i.e. by a{q = 20). The 
/(a) curves for the subregions display a variety of shapes. 

We have attempted to determine whether the subregions could be uniquely ordered according 
to some overall property of the f{a) curves. Most of these attempts consisted of plotting the 
positions of the subregions in a two-dimensional diagram, each of whose axes correspond to some 
property of the /(a) curves. For example, one such property is the range of / values covered by a 
fixed range of q values which we chose as 4 < < 20 (since the q < 4 points arc very similar for 
most regions). This parameter A/ ^ f{q = 4) — f{q = 20) may be thought of as a "dimensional 
diversity," since images dominated by a single strong point-like concentration will have / values 
close to zero for most positive q values (low dimensional diversity), while those with a range of 
geometries (from surface-like to point-like) will have a large value of this parameter. Figure 8 
shows the observed regions in the A/ — a{q = 20) plane, where a{q = 20) indicates the strength 
of the strongest concentrations. A fair correlation is present, but this is expected because the 
/(a) curve must, by its definition, remain below the diagonal line / = a, so the range in / values 
must decrease for curves which have smaller a values; i.e. images which possess strong dominant 
concentrations will usually have smaller diversities. 

The ordering of the regions in Fig. 7 does correspond to the visual appearance of the geometry 
of the structures in the column density images. For example, R CrA W has the broadest /(a) and 
has the largest "point-like" contribution, consistent with the fact that the region is dominated by 
an extremely high-column density core. L134 and Cham I are also "point-like," although not as 
concentrated as R CrA W. At the other extreme, the /(a) curves for Oph N, W, and U are the 
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narrowest of all the regions, in agreement with the smoother, diffuse, appearance and laek of strong 
concentrations seen in the column density images. The vertical extent of the /(a) curves for the 
most diffuse regions (Oph N,W,U) may be unreHable because of finite size effects. For example, 
it is well-known that a white noise field, which should appear as a single point at a = / = 2, 
will instead exhibit a very narrow measured f{a) that extends downward to smaller dimensions 
because of the finite size of the sample. Thus the dimensional diversity f(q=4)-f(q=20) of Fig. 8 
could be significantly smaller for these regions than calculated here. The intermediate regions 
(Musca, p Oph, R CrA E, Cham II, Cham III, Taurus, and Oph E) all exhibit some quasi-linear 
or filamentary structure, mixed with point-like concentrations and more diffuse structures. Oph 
E does not contain any prominent filaments, being primarily diffuse like Oph N,W,U, but it does 
contain an irregular "ridge" that may reflect the presence of an ionization/shock front, and it is 
probably this feature that shifts the /(q) spectrum toward the regions with prominent filaments. 
Notice also that the total R CrA region visually resembles the p Oph region dense core with 
"streamers" to one side. It is because we have divided R CrA into two subregions (the core and 
the streamers) that the core, R CrA W, separates so clearly from the rest of the regions in Fig. 8. 

This result illustrates the problematic dependence of all structural descriptors on the 
subjective pre-selection of the region to be described. 




- 



1 21 21 21 21 21 2 
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Fig. 7. — Spectrum of singularities of the column density maps of all the subregions. Due to 
background and calibration uncertainties only the left positive-g portion of the /(a) curves were 
calculated. The regions are ordered by the strength of the strongest singularities present, i.e. by 
a{q = 20) which is the width of the /(q) curve. 



-14- 



1,5 



o 



I 

II 

cr 



.5 



1 — 1 — r 



T [ I I [ I 1 I I |— J— |— I I |- 

Oph N 



Oph W 



Taurus 



Oph U 



pOph 



Oph E 



RCrAE 



L 134 



Cham II 



RCrAW 



Musca Cham 111 



Cham I 



..l.-„,i L 



I I I I ' I I \ I I L 



I I I 1 I L 



1.2 1.4 
a ( q = 20 ) 



1.6 



1.8 



Fig. 8. — The "dimensional diversity," taken as f{q = 4) — f{q = 20), is plotted as a function of 
the strength of the strongest singularities present in the region, i.e. a{q = 20). 

Since the ordering of regions described above does involve the degree of concentration, it 
might be erroneously thought that the properties of /(o) at large q are simply reflecting the degree 
of density contrast. We have examined the probability distributions of column density for all the 
subregions and their moments. The standard deviation (s.d.) relative to the mean column density 
varies from 0.20 (Musca) to 0.55 (R CrA W). So R CrA W is extreme in both f{d) ordering and 
column density contrast, and, similarly, Oph U and Oph N arc rather extreme at the other end, 
ranking 10 and 8 (out of 11) in standard deviation, respectively. However there are exceptions 
that make it clear that the f{a) ranking is primarily geometrical, even if density contrast does 
correlate with geometry in some cases. For example, Musca has by far the lowest s.d., yet lies in 
the middle of the regions in Fig. 9 because it is primarily a simple filament. L 134 is ranked 7 in 
s.d. (0.34), yet it, along with Cham I, is, after R CrA W, the most "point-like" of the regions, 
not at all like Oph N, which has a similar s.d. Another example is R CrA E, which is ranked 2 in 
s.d. (0.46), even though it is near Musca (rank 11), Cham II (rank 9), and Cham III (rank 6) in 
Fig. 8 due to the prevalence of line-like structure in all these regions. Oph W is among the three 
most "space-filling" or homogeneous regions according to Fig. 8, yet its standard deviation is rank 
5 (0.38). Clearly the /(a) indicator is responding to geometry and not density contrast in these 
cases. 

It should be emphasized that, while some properties of /(a) do reflect some of the gross 
character of the geometry of these regions, and hence may be useful descriptors for classification 
and comparison with simulations, f{a) definitely docs not capture other aspects of structural 
detail. For example, the Cham II and Cham III regions have very similar /(a) curves but different 
visual appearances: The Cham III region looks more "disrupted" or filamentary (even though 
it contains few, if any, young stellar objects; see below). Thus, the /(o) spectrum, like all other 
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structural descriptors that have been proposed, suffers from severe compression of structural 
information. On the other hand, the utihty of different structural descriptors, even if, or just 
because, they are highly compressed, is their ability to capture different aspects of the structure, 
and in that regard the geometrical information expressed compactly by /(a) is important because 
those aspects are not expressed by other descriptors. However it appears that no quantitative 
descriptor will be capable of capturing any but a small fraction of the aspects of structure apparent 
to the eye/brain. This is of course a common and severe problem in all applications of pattern 
analysis when description, rather than specific recognition or identification, is the goal. 

Our method of ordering the regions in terms of some property of the /(a) curve is not unique, 
and we considered other approaches. For example, we tried to order the regions according to the 
distance of /(a) from some fiducial reference /(a) curve, defining distance as in Adams (1992) and 
Adams and Wiseman (1994). For the fiducial curve we tried various colored noises, stochastic 
multiplicative processes, and the mean /(o) for the observed region. However it was found that 
the resulting ordering depended strongly on the choice of the standard /(a) curve. We suspect 
that this dependence of ordering on the adopted reference function will be a general problem 
for this approach to classification. Instead we would suggest that the variation of ordering with 
respect to fiducial reference function is itself an important descriptor of complexity, since it is an 
indication of the number of ways in which the collection of regions can be viewed; i.e. a measure of 
the ambiguity presented by the structure. This has been a common theme in recent discussions of 
the meaning of complexity, although we know of no quantitative measure that has been proposed. 

The relation of f{d) to the stellar content of the different regions is also of interest. It is 
generally thought that the star formation properties of cloud complexes might be related to the 
geometry of the dust and gas, although even the direction of the causation is unclear. Larson 
(1992) proposed that the stellar initial mass function may be a consequence of the geometry of the 
gas from which stars form. From a different point of view, it is generally thought that the presence 
of young stars should alter the structure of the gas in the parent cloud complex in a severe manner, 
although no quantitative descriptor of this effect has, to our knowledge, been proposed. For this 
reason, several star formation indicators were compiled for each of the regions, including the 
surface density of young stellar objects, the ratio of stellar luminosity to gas mass, and the star 
formation efficiency. (See the references in Table 1; also the study of star formation in these regions 
by Chen et al. 1997.) No correlation with the shape of the /(a) curve was found. The reason for 
this lack of correlation is fairly obvious: There is simply no one-to-one correspondence between 
gas morphology and stellar content. For example, while L134, with little internal activity, has a 
much smaller column density than either p Oph or R CrA W, it certainly does have a "core-like" 
appearance, and so is geometrically similar to the cores with strong star iformation activity. When 
discussing "cores" , this lack of correspondence may seem trivial, since it might simply be due to 
the fact that a region without internal activity hasn't had time to begin star formation, or has too 
small a mass to be bound, etc. However other examples show that the lack of correspondence goes 
further than regions that appear to comply with conceptions of quasi-spherical "clouds." Both 
Taurus and Cham III have similarly prominent filamentary structure, but Cham III is devoid of 
star formation activity. Musca is another filament with little internal activity; the same applies to 
the streamers extending away from the p Oph and R CrA cores. The differences may be due to the 
ages of the regions or because the gas morphology reflects external forces (particularly for Musca) 
as well as internal processes like self-gravity or hydrodynamic instabilities. Perhaps the most 
interesting comparison is Cham II and III. As mentioned above, Cham III appears significantly 
more "disrupted," even though it has essentially no measurable internal star formation activity. 
These examples suggest that it will be quite difficult to understand any relationship between star 
formation and gas structure, given the small number of regions that can be mapped at sufficient 
resolution, and the fact that we observe each region at different effective ages. Still, it would be 
interesting to carry out the type of geometrical analysis presented here for the column density 
structures in regions with strong massive star formation. The IRAS data are unsuitable for such 
maps because the heating by the massive stars probably seriously distorts the derived column 
densities, as explained in sec.II. The Bell Labs ^^CO survey provides sufficiently dense mapping 
of massive star formation regions for estimates of /(a) (for the column densities and the radial 
velocities), and should be attempted. Based on the results found here, it is likely that /(a) will 
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primarily reflect the mix of strong "point-like" structures (very dense "cores") and filamentary 
structures, both of which are known to exist in these regions. Thus, there is no reason to expect 
the properties of the f{a) curve to be significantly different from those found for, say, Taurus or p 
Oph. 



5. Summary 

We have investigated the possibility of multifractal scaling in several regions and subregions 
of local cloud complexes using column density images with a very large spatial dynamic range 
constructed from IRAS data. The major results are as follows. 

1. All the regions and subregions exhibited multifractal scaling over the range 0.4 to 4 parsecs, 
the largest range available for investigating scaling using the available images and this technique, 
even though the images covered a much larger range, because of the sensitivity to finite size and 
resolution effects (see also Takayasu and Suzuki 1991). That the column density images do exhibit 
multifractal power law scaling is a new and surprising result, which should provide an important 
constraint on theoretical models for the physical origin of the structure in star- forming clouds. In 
particular, since MFS is commonly associated with multiplicative processes, the existence of MFS 
in cloud complexes can be taken as indirect evidence for nested hierarchical structure. We have 
shown that the correlation function of the a-images of Taurus exhibits a distinct range in which 
logarithmic decay occurs, consistent with predictions for a hierarchical multifractal process. 

2. Contrary to results obtained for the dissipation fields and passive scalar fields of 
incompressible turbulence (Sreenivasan 1996 and references therein; also Muzy et al. 1991), the 
/(a) curves for these interstellar regions have a variety of shapes indicating that they do not belong 
to a single universality class with respect to f{a) . This result suggests to us that supersonic 
magneto-gravitational "turbulence" docs not possess a statistical equilibrium range, or that any 
such quasi-equilibrium evolves significantly with time or depends on initial or environmental 
conditions. This is in contrast to the similar perimeter-area dimension found in previous studies 
of local star- forming cloud complexes. However it remains to be seen whether the radial velocity 
images of these regions are multifractal, and whether they exhibit a universal MFS. 

3. The forms of the f{a) curves for these regions do seem to correlate fairly well with the 
mixture of geometric forms seen visually. Thus the MFS appears to be a useful descriptor for 
capturing this aspect of structure, and should be useful in combination with other descriptors in 
eventually obtaining a quantitative classification scheme for interstellar structure. However we 
find no correlation between the forms of the /(o) curves and the level of internal star formation 
activity. This lack of correlation between morphological structure and star formation activity can 
be easily seen directly in several examples by eye. 

The fact that all the regions do exhibit well-defined multifractal scaling provides a new 
constraint on theoretical models for cloud evolution and star formation. Simulations of realistic 
ISM evolution are approaching sufficient resolution to allow a search for multifractal scaling, 
although no such search has yet been conducted. It is known that simple models for both 
hierarchical fragmentation (e.g. sequential probabilistic partitioning, as has been used to model 
the f{a) spectrum of incompressible turbulence) and coalescence (aggregation) can result in 
multifractal scaling. However, in the fragmentation case the dominant physics is unclear, since the 
structure might be produced by nonlinear fluid advection, as in incompressible turbulence, or by 
repeated gravitational instability or other fluid instabilities. In the case of coalescence, multifractal 
scaling has only been demonstrated for Monte Carlo simulations in which spatial correlations were 
built in using a spcciflc power law correlation function (Menci ct al. 1993). It is hoped that the 
present results will motivate theoretical searches for multifractal scaling in sinuilations of cloud 
evolution based on the hydrodynamic equations. Although very large spatial dynamic range is 
required, even a negative result would allow one to conclude that an important physical process 
has been neglected, while a positive result might allow an identification of the physical processes 
that are responsible for the variation in /(o) found among different star-forming regions. 
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A. Validity of IRAS Column Densities 

In order to understand the degree to which the IRAS 100 /xm column densities, Nioo, 
represent the true column density distribution, we have compared our Nioo images to images 
of five subregions in Taurus constructed from ^^CO (Heyer et al. 1987), visual extinction Ay 
(Ccrnicharo and Bachillar 1984), and, in the case of L1495, C^'^^O (Duvert et al. 1986). These 
subregions arc L1529(=B18), B217, L1495, L1536, and Heiles Cloud 2. Both ^^CO and Ay were 

available for Hcilcs Cloud 2 and L1529, while only ^^CO was available for B217 and only Ay for 
L1506, L1495, and L1536. 

A large-scale logarithmically scaled Nioo image of Taurus is presented in Figure 9. The image 
has been processed by unsharp masking to bring out the smaller-scale structure. The unsharp 

masking is equivalent to artificially modifying the column density power spectrum by putting more 
power at larger wavenumbers. We were careful to verify that this image enhancement did not 
introduce any spurious or artifactual features. The calculations in the text and the comparisons 
described below are based on unprocessed images. 

The comparisons all yielded surprisingly positive results in the sense that the resulting 
structures are very similar in visual appearance. A comparison of ^^CO, Nioo, and Ay for the 
L1529 subregion (the filament in the lower center of Fig. 9) is presented in Figure 10. The 

^^CO and Ay images were taken in digitized form directly from the original data, but have 
been smoothed to the adopted resolution of the IRAS data. The agreement in morphology is 
good, especially compared to the ^^CO data. Notice the nearly circular ring extending from the 
middle northern part of this substructure, which can be seen to various degrees in all three maps, 
although its ring-like form is barely discernable in the extinction image. We concur with Heyer 
et al. 's conclusion that the extinction errors must be considerably larger than usually quoted in 
the literature. 

Fig. 10 is typical of the five comparisons that were carried out. The only discrepancy was 
found for the very opaque core of L1495. Luckily, a C^^O map was available for this region (Duvert 
et al. 1986), which showed that the "extra structure" seen in the IRAS map does indeed appear 
similarly in the C^^O map; the agreement is only compromised by the presence of an embedded 
source which produces a spurious "hole" in the IRAS map-otherwise the agreement is surprisingly 
good. We conclude that the ^^CO and Ay maps are saturated, and that the IRAS 100 /xm column 
density map probes optical depths similar to that probed by C^^O, at least in this case. 

Another indication that the IRAS column density maps presented here faithfully represent 
the structure is the excellent agreement between our images and the adaptive grid/ wavelet 
transform optical extinction maps presented by Cambresy (1999) for Taurus, Sco-Oph, R Cr 
A, and Cham-Musca. The techniques used by Cambresy allow extinction to probe much larger 
column densities, and more accurately, than previous extinction studies. 
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Fig. 9. — 100 fim column density map of Taurus, with unsharp masking applied to enhance small 
scale structure. 



Fig. 10. — Comparison of column density images for L1529 in Taurus based on ^^CO (top), Nioo 
(middle), and Ay (bottom). 



B. Multiplicative Toy Models 

The f{a) spectrum is ideally suited to characterizing structure which can be created by 

a multiplicative process. We present one such process, called the p-modcl (see Sreenivasan 
1991) or multinomial multiplicative process (Feder 1988), to demonstrate the sensitivity of the 
/(a) spectrum to both scaling and geometry in images. 
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Fig. 11. — The multiplicative construction process defining the toy models. A square with an 
initially uniform measure is divided into four smaller squares. In each, the original measure is 
multiplied by a constant Pj, where ?e{0, 1,2,3}. The set {Pi} remains fixed throughout the entire 
construction process. Each of the squares are then divided once again, and the construction process 
is repeated a large number of times. The final function will consist of a large number of spikes and 
is known to possess multifractal scaling. This is Sreenivisan's (1991) "p-model." The /(a) curves 
for several multifractals generated with the p-model are shown in Fig. 12. 

Figure 11 shows the construction process for the model. A uniform square is divided into 
four smaller squares each receiving a probability p\. Each smaller square, in turn, is subdivided 
into still smaller squares, the probability in each given by the product of pi with the probability 
of the "parent" square. Iterating this process produces a probability distribution densely filled 
with singularities. Singularities with common exponents, a, are distributed as homogeneous 
fractals; however the fractal dimension varies as a function of the exponent (Feder 1988). The 
/(a) spectrum provides a convenient method of characterizing both the scaling exponents of the 
singularities and the fractal dimension of the "iso-a" sets. 
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Fig. 12. — The /(a) singularity spectrum for the toy multiphcative models defined in Fig. 11. In 
the right-hand panels the diagonal probabilities are set equal P\ = P2 and Pq = P^ = 0.5 — -Pi. In 
the left-hand panels three probabilities are set equal Pi = P2 = P3 with Pq = 1 — 3Pi. In each 
case, the difference |Pi — Pq\ increases up the page producing stronger "singularities." 

Figure 12 shows the /(a) spectrum for several p-models. In each case the pattern of 

probabilities was kept constant for all stages of the construction. In the left hand column, 
the largest probability is given by pi and the remaining three probabilities are set equal 
7^2 = P3 = P4 = (1 — Pi)/3- The ratio pi/p2 is decreased down the column. All the curves peak 
at /(q' = 0) = 2 corresponding to the two-dimensionality of the images. The asymmetry of the 
/(a) curves reflects the asymmetry in the spatial distributions of the strong and weak singularities; 
the strongest singularities arc distributed as isolated points {f{q = 00) = 0), whereas the weakest 
singularities are more densely distributed and can be shown to lie on a "Serpinski gasket" with a 
fractal dimension f{q = —00) ^ 1.58. 

The value of a at the q = 00 point represents the strength of the strongest singularities and 
decreases as the ratio pi/p2 increases. The width of the /(a) curve measures the range of the 
singularity strengths present in the image. 

The curves in the right-hand column are based on p-modcls in which pi = p^ and 
P2 = Pa = (1 ~Pi)/2- Again, the ratio pi/p2 is decreased down the column. The strongest 
singularities occur where the probability distribution is largest, i.e. where the probabilities pi or 
P3 have been used in each stage of the construction process. Since the number of such squares at 
the n^^ construction stage scales as n oc L~^, where L is the box size, the strongest singularities 
have a fractal dimension D = 1. Similarly, the weakest singularities also have D = 1, producing a 
symmetrical f{d) curve which remains above / = 1. 

We conclude with a cautionary remark. Figure 13 shows three point sets which evoke different 
visual impressions yet which possess the same fractal dimension of D = 1.0. It is obvioTis that 
the fractal dimension (and therefore the f(a) spectrum) is a severely incomplete and degenerate 
description of the underlying structure and care should be taken when making inferences based on 
this (or any other) structural descriptor. A well-known similar example is the power spectrum 
(and hence correlation function), which is highly degenerate (very dissimilar structures can 
produce nearly-identical power spectra) because phase information is lost in its construction. 
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More detailed discussion of the ambiguity involved in interpreting the multifraetal spectrum can 
be found in Chhabra, Jensen, and Sreenivasan (1989) and Bulajich and Perez-Pascual (1991). 
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Fig. 13. — Examples showing the degeneracy of the fractal dimension as a structure descriptor. 
Each point set has a fractal dimension D = I. The fractal in the center panel was created with the 
fixed p-model with Pi = P2 = 0.5 and Pq = P3 = (see Fig. 11). The ordering of the probabilities 
was reversed on alternate stages of the construction process. The fractal in the right-hand panel 
was constructed using the same set of probabilities, but with random ordering at each stage of the 
construction process. 
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Table 1. Properties of the Regions Studied 



Subregion 


R. A. 


Dec. 


# Res. Elem. 


Dist. (pc) 


# YS0/pc2 


Tau 


4'^30' 


+26°.0 


3.2 X 10^ 


140'' 


0.3^ 


RCrA W 


19'*04' 


-37°. 


1.0 X 10^ 


130^ 


3.0*^ 


R CrA E 


19''32' 


-37°.0 


2.6 X 10"^ 


130'= 


0.0*^ 


Cham 1 


ll'^OO' 


-77° .0 


1.7 X 10"^ 


140*= 


3.5/ 


Cham 2 


12'*40' 


-78°.0 


1.1 X 10^ 


140^ 


2.0^ 


Cham 3 


12^^20' 


-80°.0 


1.4 X 10"^ 


140^ 


0.0-^ 


Musca 
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140^= 


O.l-'' 


p Oph 
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-26°. 


5.6 X 10^ 


125f 


23.0'' 
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16'*50' 


-24°. 


2.0 X 10^ 


1255 




Oph W 


15''40' 


-27°. 


3.6 X 10^ 


125^? 




OphU 


16^26' 


-19°.0 


2.2 X 10^ 


125^^ 




OphN 


16'*40' 


-10°.0 


1.0 X 10^ 


1259 




L 134 




-5°.0 


2.4 X 10^ 


160^ 


0^' 
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